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Abstract 

The short-range and tensor correlations associated to reahstic nucleon-nucleon interactions in- 
duce a population of high-momentum components in the many-body nuclear wave function. We 
study the impact of such high-momentum components on bulk observables associated to isospin 
asymmetric matter. The kinetic part of the symmetry energy is strongly reduced by correlations 
when compared to the non-interacting case. The origin of this behavior is elucidated using realistic 
interactions with different short-range and tensor structures. 

PACS numbers: 21.65.Ef: 26.60.-c 
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The existence of high-momentum components in the nuclear many-body wave function 
is a well-estabhshed property both from an experimental and a theoretical points of 

view I2, 4|_. Short-range correlated pairs have been studied in detail at electron scattering 
facilities [5]. Two-nucleon knock-out reactions have identified the predominance of isospin 
/ = correlated pairs [2], an effect that has been related to the tensor component of the 
nucleon-nucleon (NN) interaction js, 7|. 

The enhanced effect of correlations on neutron-proton (np) pairs with respect to neutron- 
neutron (nn) pairs suggests that correlations can be increased (or even tuned) in isospin 
asymmetric systems, where a different number of np and nn pairs exist. Our microscopic 
many-body calculations indicate that the dependence of short-range and tensor correlations 
on the isospin asymmetry of the system, a = [N — Z)/{N + Z), is well understood from 
general theoretical principles jsj. One-body occupations, for instance, follow a systematic 
trend: neutrons become less depleted as the system becomes more neutron-rich, while the 
proton depletion increases js, 9|. 

Here, rather than looking at microscopic properties, we want to quantify the effect that 
NN correlations have on the bulk properties of isospin asymmetric systems. By analysing the 
energy of symmetric, asymmetric and neutron matter obtained within a realistic many-body 
approach, we will draw conclusions on the impact of correlations on asymmetric systems. 
We will focus our attention on the symmetry energy, which characterises the properties of 



isospin-rich nuclei as well as neutron stars 
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To quantify the effect of correlations in a meaningful way, we use a many-body approxima- 
tion that includes consistently short-range and tensor correlations. The ladder approxima- 
tion, implemented within the self-consistent Green's functions (SCGF) approach, provides a 
microscopic description of these effects via a fully dressed propagation of nucleons in nuclear 
matter ^. This is achieved by (a) computing the scattering of particles via a T- matrix (or 
effective interaction) in the medium, (b) extracting a self-energy out of the effective interac- 
tion and (c) using Dyson's equation to build two-body propagators which are subsequently 
inserted in the scattering equation [sl, To solve the close set of equations, an iterative 
numerical procedure is a must. Recent advances have allowed implementations both at zero 
[1^ and at finite temperature [isl using fully realistic NN interactions. We will focus here 
in calculations based on two-body forces only, including partial waves up to J = 4 (J = 8) 
for the dispersive (Hartree-Fock) contributions. 



2 



Fully self-consistent ladder calculations in isospin asymmetric nuclear matter are scarce 
[sl, qI- Previous results have mostly highlighted the different role of correlation in either 
symmetric or pure neutron matter [11]. We will analyse here the asymmetry dependence of 
the different components of the energy. More specifically, we want to look into the kinetic 



component of the energy. The latter has been identified recently as a particularly good indi- 
cator of correlations in the nuclear ground state [l4, [isj. We will confirm these indications 
using SCGF techniques. 

The bulk properties of nuclear and neutron matter are obtained within the SCGF ap- 
proach via the Galitskii-Migdal-Koltun sum-rule: 



where i/ = 4 (2) is the degeneracy of nuclear (neutron) matter, p is the total density and 
/(w) = [1 + exp (w — ;u)/T]~^ is a Fermi-Dirac distribution. Roughly speaking, the one- 
body spectral function, A{k^uj\ represents the probability of knocking out or adding a 
particle with a given single-particle momentum, fc, and energy, w. The spectral function also 



gives access to all the one-body operators of the system 16|. For instance, the momentum 



distribution, n[k\ is obtained by convoluting the spectral function with a Fermi-Dirac factor: 

n{k) = I ^A{k,co)f{u). (2) 

Correlations beyond the mean-field approximation have a particularly clear manifestation 
in the momentum distribution [9]. A sizeable dep letion appears below the Fermi sea, while 



high-momentum components are populated 



171]. To illustrate this point, we show in 



the upper left (right) panel of Fig. [T] the momentum distribution for symmetric nuclear 
(pure neutron) matter. The results obtained within the SCGF method for the Argonne vl8 



(Avl8) 18| and the CDBonn 19| interactions (solid and dash-dotted lines, respectively) are 
compared to the momentum distribution of the Free Fermi Gas (FFG) in the same condi- 
tions (dashed lines). The FFG is used here as a benchmark for thermal effects. A common 
feature of the SCGF and the FFG results is the softening of the distribution around the 
Fermi surface, k = kp, associated to the finite temperature of T = 5 MeV. Calculations are 



performed at this temperature to avoid the pairing instability 20|. The density, p = 0.16 
fm~^ does not correspond to the saturation point of any of the two potentials. For Avl8 
(CDBonn) within the SCGF approach at T = 5 MeV, the saturation density is at p = 0.19 
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FIG. 1. (Color online) Upper left (right) panels: momentum distribution of symmetric nuclear 
(neutron) matter obtained with the SCGF approximation for Avl8 (full lines) and CDBonn (dash- 
dotted lined) and for the FFG (dashed lines). Lower panels: the momentum distribution times 
k'^/kp, a quantity directly connected to the kinetic energy. All the results are computed at p = 0.16 
fm ~^ and T = 5 MeV. The dotted line shows the high momentum constant of Ref. [2^ 



(0.32) fm . Three-body forces (3BFs) should improve the saturation properties without 



strongly affecting isovector properties 12j. 



Correlation effects in the momentum distribution are substantially different in symmetric 
nuclear matter (SNM) and in pure neutron matter (PNM) j^. The presence of the S — D 
tensor interaction in the isospin-saturated system induces a larger amount of correlations. 
As a consequence, the Fermi surface is quite more depleted for SNM than for PNM (compare 
the upper left and right panels in Fig. [1]). Typical values for these depletions are obtained 
from the occupation at zero momentum: n(0) ~ 0.87 for SNM and n(0) ~ 0.96 for PNM. 
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CDB 


kp 
kp 2kp 
2kp oo 

oo 


0.762 15.8 -12.6 
0.211 12.0 -1.54 
0.027 6.95 -0.59 
1.00 34.7 -14.7 


0.869 28.9 10.1 
0.121 9.87 3.78 
0.010 3.78 0.28 
1.00 42.5 14.1 


Avl8 


kp 
kp 2kp 
2kp oo 

oo 


0.755 15.6 -7.65 
0.194 11.4 -0.997 
0.051 14.5 -1.29 
1.00 41.5 -9.94 


0.863 28.7 11.6 
0.119 10.3 3.24 
0.018 7.16 0.32 
1.00 46.2 15.2 


FFG 


kp 
kp 2kp 
2kp oo 

oo 


0.861 17.7 17.7 
0.139 6.00 6.00 
0.00 0.00 0.00 
1.00 23.7 23.7 


0.912 30.4 30.4 
0.089 5.75 5.75 
0.00 0.00 0.00 
1.00 36.2 36.2 



TABLE I. Contributions of different momentum regions to the total density (columns 3 and 6), 
kinetic (4 and 7, in MeV) and total energies (5 and 8, in MeV) for SNM (columns 3, 4 and 5) and 
PNM (columns 6,7 and 8) with different NN interactions. The FFG case is also included. All the 
results are computed at p = 0.16 fm^^ and T = 5 MeV. 

Because the momentum distribution is normalised to the total density, the high momen- 
tum components are also rather different for both systems at the same density. One can 
summarise these differences by looking at the integrated strength over different regions of 
momenta: 

<Pm{h, kf) = ^ r dA; k^n{k) . (3) 

The integral with m = 2 represents the fractional contribution of a given momentum region 
to the total density. Similarly, the integral with m = 4 is related to the total kinetic energy 
of the system. 

In Table [H we give the integrated strengths of SNM and PNM at p = 0.16 fm"^ and T = 5 
MeV for the Avl8 and the CDBonn interactions. As expected, in SNM there is a substantial 
depletion of states below the Fermi surface, i.e. only ~ 75 % of the strength is in the region 
between and kp. Part of the depletion has a thermal origin, and the comparison with the 
FFG in the same momentum region suggests that between 1/2 and 2/3 of the integrated 
depletion comes from the softening of the Fermi surface due to the finite temperature. The 
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effect of correlations is also important beyond the Fermi surface: for SNM (PNM) there is 
still a 3 — 5% (1 — 2%) of the strength in the region k > 2kF- The energy per particle is 



also affected by short-range correlations 21|. In particular, the kinetic energy, 

-n{k) , (4) 



A p J (27r)3 2m 

increases with respect to the FFG due to the population of high-momentum components 
We present the contributions of the different momentum regions to the kinetic energy 
of nuclear and neutron matter in columns 4 and 7 of Table [H First of all, let us note that the 
total integrated values of the correlated kinetic energies are larger than those of the FFG. 
For SNM, these are ~ 11 and ~ 17 MeV larger for CDBonn and Avl8, respectively. The 
difference between the two potentials is expected, since Avl8 is a hard interaction and has a 
stronger tensor component than CDBonn. In PNM the total kinetic energy is 6 and 10 MeV 
higher, respectively, for both interactions. This emphasises again the idea that correlation 
effects play a smaller role in PNM than in SNM. 

Let us stress once more the importance of components beyond kp in the kinetic energy. 
For SNM, these amount to more than 50 % of the total, while in PNM they account for more 
than 25 %. In contrast, for the FFG at this temperature, the contribution of states above 
kp is less than 25% (15%) for SNM (PNM). The FFG strength above the Fermi surface is 
due to thermal effects, which tend to be localised within a small region around kp for low 
temperatures. As a consequence, there are no contributions to the FFG energy in the region 
oik > 2k p. The contributions in this region in the interacting case can be entirely attributed 
to NN correlation effects. For a hard interaction like Avl8, the contribution beyond 2kF in 
SNM is even larger than that between the Fermi surface and 2k p. 

A visual representation of the contributions of different momentum regions to the kinetic 
energy is obtained by looking directly at the integrand in the formula for the kinetic energy, 
k'^n{k). The lower panels of Fig. [1] show this quantity for SNM (left) and PNM (right). 
Compared to the FFG, one finds that the integrands for SNM and PNM have substantial 
contributions for k > kp. In both cases, the integrand extends to very high momenta, up to 
4 — 5kp. With the present normalization, the high momentum components of both systems 
are similar in size, but they extend to higher momenta for SNM than for PNM. As already 
discussed, the kinetic energy of SNM is higher than that of PNM with respect to the FFG. 
This will have a strong impact on the kinetic component of the symmetry energy. 
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If SNM or PNM were in the unitary regime, the integrands shown in the lower panels of 
Fig. [1] would tend to the so-called contact constant, C/lNkp), at very high momenta 
For a unitary gas, the momentum distribution in the /c — ?■ oo region would therefore decay 
as k~^. SNM and PNM are not unitary gases, though, and their kinetic energy integrands 
are not constant at high momentum. Instead, these functions have a peak at kp and then 
decreases sharply (at T = there would be a discontinuity at kp). In the region 1.5 — Skp, 
k'^n{k) levels off (signalling a plausible n{k) ~ scaling in this region) and then decays 
softly as k increases. The presence of a secondary maximum in the Avl8 results explains 
why the kinetic energy has such a large contribution in the k > 2kp region. To ensure 
the correct high momentum limit, the integrands have been normalised by the constant 
(37r^)/(z/), i.e. the normalizations in SNM and PNM differ by a factor of 2. With this, 
the magnitude and shape of the high momentum components look quite similar for the two 
systems and interactions considered in Fig. [H For SNM (PNM) under these conditions, 
we have T/sp ~ 0.14(0.08) and we find that, in both cases, the maximum value of the 
integrand is below 2. In the same conditions, quantum Monte-Carlo calculations of a unitary 



gas suggest a value of C/{Nkp) ~ 3 [22 1. 



To compute the symmetry energy (and its kinetic and potential components), one gener- 
ally resorts to the parabolic formula, i.e. one assumes that the energy per particle (or any 
of its components) has a quadratic dependence on asymmetry, 

j{p,a) = j{p,0) + S{p)a\ (5) 

This immediately yields that the symmetry energy, S{p), is given by the difference of PNM 
and SNM energies: 

S{p) = j{pA)-j{p,0). (6) 

The SCGF approach can be generalised to isospin asymmetric systems [8|. We can thus 
validate the parabolic assumption by performing calculations at different a's and looking at 
the explicit dependence of the energy on isospin asymmetry. We have performed this check 
for the three components of the energy (kinetic, potential and total) and two representative 
potentials, CDBonn (circles) and AvlS (squares), at p = 0.16 fm"'^ and T = 5 MeV. The 
results are displayed in the three panels of Fig. |2l as a function of a^. Note that the vertical 
extent of all three panels is the same. In general, the three components seem to have a 
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FIG. 2. (Color online) Isospin asymmetry dependence of different components of the energy for the 
CDBonn (circles) and Avl8 (squares) potentials. The upper, central and lower panels correspond 
to the kinetic, total and potential energies, respectively. The triangles of the upper panel give the 
energy of the Free Fermi Gas in the same conditions, p = 0.16 fm~^ and T = 5 MeV. Dotted lines 
are linear regressions to guide the eye. 



well-defined parabolic dependence on a. To stress this point, we show the results of a linear 
regression fit (dotted lines) on top of the calculated data. The slope of the linear regression 
reduces to the different components (kinetic, potential and total) of the symmetry energy if 
the parabolic approximation holds exactly. We have found a very good agreement (generally 
within less than 0.5 MeV) between the slopes and the values obtained using the differences 
of Eq. ^. 
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In the upper panel of Fig. [2l we compare the kinetic energy of the two potentials to that 
of the corresponding FFG (triangles). As expected, the correlated kinetic energy is larger 
than the FFG at all asymmetries. A hard interaction (Avl8) leads to a substantially larger 
kinetic energy than a soft one (CDBonn). Moreover, the isospin dependence of both forces 
is different (and different than the FFG). While the kinetic energy of Avl8 in SNM {a = 0) 
is K/A ~ 42 MeV and that of PNM is ~ 46 MeV, for CDBonn these two quantities are 
35 and 43 MeV, respectively. In other words, the difference between the kinetic energies 
of PNM and SNM is smaller for Avl8 than for CDBonn. Correspondingly, both of these 
differences are smaller than that associated to the FFG. 

The small value of the kinetic symmetry energy in a correlated approach is one of the 
major conclusions of this paper. One can understand the origin of this behaviour from the 
following reasoning. The tensor component of the NN force, acting on SNM, induces large 
correlations and produces an important renormalization of the kinetic energy with respect 
to the FFG. The absence of this component in PNM reduces the relative enhancement of 
the kinetic energy. Consequently, the difference in total kinetic energies is smaller for the 
correlated case than for the FFG value. In turn, this implies that the kinetic symmetry 
energy is reduced. Within this picture, tensor correlations, and the renormalization they 
induce in the kinetic energy, seem to be the main responsible for the small values of the 
kinetic symmetry energy. 

The asymmetry dependence of the total energy per particle is driven by a competition 
between the kinetic and the potential terms. The size of both contributions is density de- 
pendent but, at p = 0.16 fm~^, the potential term largely dominates the isospin dependence. 
This can be directly seen in Fig. [2) the difference between the PNM and SNM potential en- 
ergies is of the order of 20 MeV, while for the kinetic term the differences are below 10 MeV. 

The values of the different components of the symmetry energy obtained in the parabolic 
approximation at T = 5 MeV are given in Table [Tll In addition to AvlS and CDBonn, we 
consider two other representative NN interactions: the somewhat hard Nijmegen 1 (23! and 
the very soft N3L0 with a A = 500 MeV cut-off 2^]. Note that the calculations have been 
performed at p = 0.16 fm"'^, which is not the saturation density of any of these potentials. 
Based on the comparison with the FFG (see Fig. [3]), the effect of finite temperature in the 
symmetry energy should be rather small. At the density we are considering, the symmetry 
energy of the FFG increases from 12.4 MeV at zero temperature to 13 MeV at T = 5 MeV 
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Slot Skin Spot L 



Avl8 25.1 4.9 20.2 37.7 
Nijl 27.4 4.6 22.8 48.5 
CDBonn 28.8 7.9 20.9 52.6 
N3L0 29.7 7.2 22.4 55.2 

TABLE II. Total, kinetic and potential contributions in MeV to the symmetry energy at p = 0.16 
fm"^ and T = 5 MeV for different NN interactions. The last column gives the L coefficient in MeV, 
related to the density dependence of S'tot in the same conditions. 



(see below for further discussion). 

The total symmetry energies predicted by the SCGF approach range between 25 and 30 
MeV, just below the currently accepted ~ 32 MeV value [25|. The inclusion of 3BFs would 
bring the SCGF results closer to experiment. Although we cannot perform calculations with 
3BFs at the moment, we can estimate their importance from existing Brueckner-Hartree- 



Fock (BHF) calculations 26|. Around p = 0.16 fm ^, 3BFs tend to increase the symmetry 



energy by 3 — 4 MeV 26| . 



The symmetry energies provided by SCG 



calculations tend to be smaller than the BHF 



ones with the same two-body NN force 26|. The origin of this difference can be under- 



stood as follows. In general, the propagation of holes and the dressing of the intermediate 
propagators in the ladder equation has an overall repulsive effect in the total energy of the 
system with respect to the BHF values [l^]. This repulsive effect is larger for SNM, where 
correlations are more substantial, than for PNM. The difference between energies is therefore 
reduced and the SCGF symmetry energy becomes smaller than the BHF one. Since this 
repulsive effect increases with density, the slope of the symmetry energy as a function of the 
density is also expected to decrease. At this point, it is worth mentioning that the BHF 
symmetry energy obtained with the Avl8 interaction, at the same density and temperature 
as in Table [ni is S'tot = 28.4 MeV, i.e. about ~ 3 MeV higher than the SCGF result. Similar 
conclusions have been discussed previously in the context of SCGF calculations [111 and of 

n 

extensions of BHF theory inspired by Green's functions [27|. 

As already explained, the kinetic symmetry energy (column 3 of Table [IT!) is small. In 
most cases it lies between 4 — 8 MeV, i.e. well below the corresponding FFG value of 
S'kin = 12 MeV. Such a small value is compatible with recent observations obtained either 
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with a simplified phenomenological model, as in Ref. [l^, or with a more sophisticated 
many-body method based on the Brueckner-Hartree-Fock (BHF) theory, as in Ref. 15 |. 
Similarly, Fermi- Hypernetted-Chain (FHNC) calculations with realistic two and three-body 



interactions also yield relatively small kinetic symmetry energies |28l. |29| . One can therefore 
conclude that this effect is caused by NN correlations in the many-body wave-function. 

For illustrative purposes, we also show, in the last column of Table Ull the value of the L 
parameter, L = 3pdS/dp, obtained in the same conditions. Most SCGF results fall below 
the currently preferred value of L ~ 50 — 60 MeV 25|]. The important corrections induced 



by 3BFs will improve this result 
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26| . Compared to BHF results with the same inter- 
action, we find a decrease of the slope of the symmetry energy around p = 0.16 fm-. m 

qualitative accordance to Ref. [llj. As an example, for Avl8 at the same temperature, the 
BHF approach provides L = 47.5 MeV. 

Our results show an almost linear correlation between the values of the symmetry energy, 
S, and the slope parameter, L. Lower symmetry energies, generally associated to harder 
interactions, correspond to lower slopes. Similarly, softer interactions seem to induce higher 
symmetry energies and slopes. This confirms the trends observed in phenomenological mean- 



field calculations from a purely microscopic perspective [26|, |30|]. In other words, one can 
think of the correlation between S and L as a general property of isospin asymmetric sys- 
tems, rather than as an artifact due to the fitting procedures of mean-field parametrizations. 
We illustrate the density dependence of our results in Fig.[3l where we show the total, kinetic 
and potential components of the symmetry energy obtained with CDBonn (left panel) and 
Avl8 (right panel) at T = 5 MeV. The symmetry energy and its components grow steadily 
with density in the density range considered here. The potential component is always larger 
in absolute value than the kinetic one, thus dominating the contribution to S. The difference 
in the density dependence of the kinetic and potential contributions is subtle. While S'kin 
grows almost linearly with p, Spot seems to have a milder density dependence. In absolute 
terms, however, both contributions have a similar importance for L around saturation den- 
sity. Similar conclusions hold for other NN interactions, not shown here for simplicity. 

It is interesting to note that the kinetic symmetry energy becomes negative below a den- 
sity of about 0.04 — 0.08 fm~^. At such low densities, the effect of thermal correlations is 
expected to be important, so one might be tempted to attribute this ant i- intuitive behaviour 
to finite temperature correlations. To pin down the importance of thermal effects, we also 
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CO 0.08 0.16 0.24 0.32 0.08 0.16 0.24 0.32 ^ 
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Density, p [fm ] Density, p [fm ] 

FIG. 3. (Color online) Components of the symmetry energy for the CDBonn (left panel) and Avl8 
(right panel) potentials at T=5 MeV. Circles, squares and triangles represent the total, kinetic 
and potential contributions, respectively. The continuous (dashed) lines correspond to the FFG 
symmetry energy at T=5 (T=0) MeV. 



show in both panels of Fig. [3l the symmetry energy of the FFG at T = (dashed line) and 
T = 5 (solid) MeV. The differences are extremely small: the symmetry energy decreases 
by less than 1 MeV when going from zero to a temperature of 5 MeV in the whole density 
range. The small contribution of temperature on the symmetry energy is caused by the 
relatively similar thermal corrections of SNM and PNM jlSj. When taking the difference 
of both energies, one eliminates practically the temperature dependence. As a matter of 
fact, in both the classical (low density) and degenerate (high density) limits, the thermal 
effects on the symmetry energy disappear exactly. For the SCGF results, this suggests 
that negative kinetic symmetry energies at low densities are in fact not a thermal, but a 
correlation-dominated effect. 

A negative kinetic symmetry energy has already been found (although not explicitly men- 



tioned) in Ref. iJ], where a simplified model for high-momentum components was proposed 
by assuming a given shape of nilz). The model depends on a single, density- independent 
parameter, a, which accounts for the depletion of states below kp. The original values of 
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a seem to be somewhat extreme when compared with depletions obtained in reahstic ap- 
proaches. As a matter of fact, tuning a to reproduce our momentum distributions ehminates 
the negative values of S'kin- With zero-temperature FHNC calculations, one also obtains neg- 
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ative kinetic symmetry energies at low densities 

We would like to stress that, in principle, negative values of S'kin, unlike negative values of 
S'tot, are not associated to a thermodynamical instability [Sll- In general, our results suggest 
that one needs to take into account high momentum components to get realistic values of the 
kinetic symmetry energy. Whether or not negative values of S'kin (or their underlying cause, 
high momentum components) have an impact on either neutron star physics or transport 
simulations remains to be seen. 

We have focused here on the kinetic component of the symmetry energy because of its 
close and transparent relation to short-range and tensor correlations. The potential energy 
contribution is equally important for isospin physics. However, due to the presence of the 
mean-field contribution, high momentum correlations will be harder to disentangle in the 
potential symmetry energy. We expect both its value at saturation as well as its density 
dependence to be different than those predicted by mean-field theories. 

In summary, we have studied the isospin asymmetry dependence of the bulk properties of 
nuclear matter within the SCGF approach. We have confirmed the quadratic dependences 
on isospin asymmetry of the total, potential and kinetic energy. We have highlighted the 
effect of NN correlations and, in particular, of high momentum components by looking at the 
population of strength above kp using realistic momentum distributions. Similarly, we have 
quantified the contribution of momenta beyond the Fermi surface in the kinetic and total 
energies with the help of the Galistkii-Migdal-Koltum sum-rule based on spectral functions 
obtained in the ladder approximation. The change in nature of high momentum compo- 
nents as the isospin asymmetry is modified leads to a substantial decrease of the kinetic 
component of the symmetry energy with respect to the FFG. The tensor components of the 
NN interaction are largely responsible for this effect. The results discussed here only include 
two-body forces, but we do not expect that the overall features will change much when 3BFs 
are included. For quantitative predictions in nuclei and neutron-star matter, however, the 
effect of three-nucleon interactions is clearly needed. Work is underway in this direction. 
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